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ABSTRACT 

In this work, we extensively investigate the terrestrial planetary formation for the 
inclined planetary systems (considering the OGLE-2006-BLG-109L system as proto- 
type) in the late stage. In the simulations, we show that the occurrence of terrestrial 
planets is quite common, in the final assembly stage. Moreover, we find that 40% of 
the runs finally occupy one planet in the habitable zone (HZ). On the other hand, the 
numerical results also indicate that the inner region of the planetesimal disk, ranging 
from ~ 0.1 to 0.3 AU, plays an important role in building up terrestrial planets. By ex- 
amining all simulations, we note that the survivals are located either between 0.1^1.0 
AU or beyond 7 AU, or at the 1:1 mean motion resonance of OGLE-2006-BLG-109Lb 
at ^2.20 AU. The outcomes suggest that it may exist moderate possibility for the 
inclined systems to harbor terrestrial planets, even planets in the HZs. 
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1 INTRODUCTION 

To date, over 500 exoplanets have been discovered 
(http: / /exoplanet.eu), revealing a diversity of planetary sys- 
tems. One fascinating fact that has revealed so far is the 
population of "Hot Jupiters" - gas giants moving in very 
small orbits (periods < 8 days) about their parent stars, 
of whic h the prototype was th e first exoplanet discovered, 
51 Peg l|Mavor fc Quelodll995l ). Recently, improvement of 
measurement precisions for both Doppler and transit dis- 
covery techniques have made the detection of more than a 
dozen exoplanets in the mass range from ~ 2 to 15 
(super-Earths), among which there are also some sh ort- 
period hot super-Earths - Gl 876 d dRivera et alll2005j) . Gl 
581 e (|Mavor et al.ll2009h. CoRoT-7b (iLeger et al.l2009h and 
the Kepler- 11 system ( Lissauer et aLlhoill ). et~ These ex- 
citing observations greatly motivate us to explore the open 
issues, to understand the formation and evolution of plane- 
tary system. 

According to core-accretion model, it is generally be- 
lieved that planetary formation could be divided into sev- 
eral stages. The starting planetary formation scenario is the 
agglomeration of /jm dust g rains to produce kilometer-size d 
bodies called planetesimals (ISafronov| [l969: Wcthcrill 1980), 
however, the b uildup of part i cles from cm to me ter size 
is still unclear (|Lissaueri fl99^ : lYoudin fc Shul 12002ft . After 
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they emerge, numerous planetesimals can further form 1000- 
kilometer-sized planetary embryos (protoplanet) and even 
the cores of gas-giants via direct collision-merger process. 
In this stage, a few larger bodies would grow rapidly con- 
sidering their do minant role of gravitational focusing in th e 
runaway growth (|Safronovll 19691 : IWetherill fc Stewarull989l ). 
At the time that they are separated into various feeding 
zones, the largest bodies become massive enough to stir up 
the nearby small objects, and subsequently the oligarchic 
growth takes over. In the oligarchic stage, the largest bod- 
ies accrete surrounding masses in the disk much slower than 
they do in runaway gr owth, but still grow more quickly than 
those minor bodies (|Kokubo fc Idal 1 19981 : iGoldreich et alj 
2004). However, the oligarchic growth may cease once these 
massive bodies accrete to have their is olation masses, where 
M lso oc Af, _1/2 Ep /2 a 3 (|Armitagell2007l k and M» is the mass 
of the central star, E p , the surface density, and a, the semi- 
major axis. This leads to the formation of larger solid cores 
for jovian planets at distant orbits. Once they reach critical 
mass (~ IOMq), the solid cores begin to a ccrete gas very 
fast, and then to gr ow into giant planets (|Kokubo fc Idal 
l2002l : llda fc Linll2004t ). All aforementioned scenarios are now 
considered to be completed within 5—10 Myr, which is crit- 
ical for the b uildup of gas-giants, as inferred lifetimes of 
gaseous disks (|Haisch et al.ll200ll : IWvattll200ct ). 

The following stage is featured by terrestrial planet for- 
mation under the interference of gas-giants, which is con- 
sidered to be the longest evolution history with collisional 
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coagulation. In this late stage, the formed gas-giants, accom- 
panied with a large number of planetesimals and planetary 
embryos, have ceased migration since most of gas in the disk 
disappears. These low-mass objects, under the perturbation 
of jovian planets as well as their mutual gravitational inter- 
action, may experience a turbulent evolution process that 
may last for hundreds of million years. During this stage, 
the embryos are greatly stirred to ha ve crossing orbits, which 
bring about numerous giant impacts (jChambers fe Wetherilll 
1998). Finally, small bodies are left over after the former 
stages h ave been cleaned up and several terrestrial planets 
emer ge (|Chambersll200ll : iRavmond et alJliool IZhang fcjH 
l2009h in the system. 

Under some circumstances, the formed planets could 
be evolved into a highly-inclined configuration, rather than 
a planar disk perpendicular to the stellar spin axis, which 
seems quite natural since both the star and planets obtain 
their angular momenta from the disk. One mechanism that 
could increase the plane tary inclination i s the scenario of 
planet-planet scattering. IChatteriee et' al. (2008) show that 
after strong planet-planet gravitational scattering, the plan- 
etary inclination, in the initial range from 0° ~ 4° , could be 
dramatically excited. In their simulations, more than half of 
the final configurations have a relative inclination exceeding 
20° and nearly a quarter above 40°. Considering the tidal 
damping timescale for inclination is usually much larger 
than the star's age (|Winn et al.ll2005l ), such highly -inclined 
planets could be found in the planetary systems which have 
undergone severe gravitational scattering in their lifetimes. 

An alternate scenario for the origin of highly-inclined 
planets is that they exist in binary systems. The presence 
of a companion star can affect the orbits of planets and the 
stability of the planetary system. In tight binaries, a close 
encounter between a plane t and the companion can e ject the 
planet out of the system l|Holman fc Wiegertlll999l ). While 
in wider binaries, secular effects like Kozai resonance become 
more important. The Kozai mechanism was first introduced 
to analyze how th e orbits of i nclined asteroids were influ- 
enced by Jupiter (Kozai 1962). Now it is widely accepted 
to explain the formation mechanism of a relative highly- 
inclined configuration, such as a Jupiter- size exoplanet in 
wide inclined binaries l|lnnanen et all fl997l ^ . etc. The Kozai 
mechanism is only notable when the inclination between the 
orbits of the planet and the companion, io, is greater than a 
critica l inclination i c of around 40 ° (depending on the sys- 
tem) jThomas fc Morbidellil 119961 ). If this is the case, the 
inclination of the planet will oscillate, between i c and io, 
and is coupled to eccentricit y e through the Kozai constant 
|Thomas fc Morbidellil Il996l ) 



Hk = \/ a{l — e 2 )cosi (1) 

Takcda et al.l (|2008l ) showed that the interplay between sec- 
ular perturbations from a stellar companion and the gravi- 
tational coupling planets could trigger the chaotic growth of 
mutual inclination angles between initially coplanar planets. 

Observations have revealed that some planetary sys- 
tems have large obliquities. Recently, direct measurements of 
the planets v Andromedae (v And) c and d unveil a 30° mu - 
tual inclination ang le between them (|McArthur et alj|2010l ) . 
iBarnes et al.l (|2010h showed that planet-planet scattering is 
a plausible mechanism to explain the observed orbits of v 
And c and d, though it is unsure whether the scattering 



was caused by instabilities among the planets or by per- 
turbations from a nearby low-mass star v And B. In recent 
years, transit observations have measured A sin i for a couple 
of systems, where A is the angle between the stellar rota- 
tion axis and the orbital angular momentum of a transiting 
planet, indicatin g that some system s with large A of 30° ± 
21° for TrES-lb iiNarita et al.1120071 ). and 62° ± 25° for HD 
17156b |Narita et al.ll200Sft . Such planetary systems might 
have a chaotic dynamic al history that produce the present 
highly-inclined planets (|Chatteriee et al]|2008l ). Currently, 
it is hard to determine the proportion of the exoplanet pop- 
ulation with a significant inclination, which requires precise 
measurements of A for a large quantity of the systems. How- 
ever, in this sense, it is extremely significant for one to make 
a deeper understanding of complicated building process for 
planetary formation, by investigating the configuration of 
highly-inclined planets. 

In this work, we extensively perform a series of groups 
of simulations, to explore the planetary formation in the 
late stage for the inclined system (application to OGLE- 
2006-BLG-109L, see Table[T]for details), by considering ma- 
jor factors that might have influence on final assembly of 
terrestrial planets. The paper is structured as follows: §2 
describes numerical setup for investigation. §3 presents the 
results and the discussion is given in section §4. Finally, we 
briefly summarize the main outcomes in §5. 



2 NUMERICAL SETUP 

In our simulations, the initial conditions emulate the en- 
vironment at the beginni ng of the chaot i c phase of ter- 
restrial pla net formation ^Chambers! l200ll ; IRavmond et alj 
2004, 2006), Moon-sized planetesimals and Mars-sized plan- 
etary embryos are distributed in the disk, in company with 
the OGLE-2006-BLG-109Lb,c analogues. The orbital ele- 
ment s are adopted from the observation data (|Gaudi et al.l 
2008). For each run, we assume the initial objects to fol- 
low a distribution with orbital radius as N oc r" 1 ^ 2 , corre- 
sponding to the annul ar mass in a disk with surface density 
ilWeidenschillindl 19771 ). 



E oc r 



-3/2 



(2) 



The Hill radius Rh within which the gravity of planet dom- 
inates is, 



Rh = (M P /3M, 



(3) 



where M p and Mq are the masses of planet and central star 
respectively. Then the feeding zone of a planet at a given 
position a can be expressed as, 



F a = CR H 



(4) 



where C is a constant. As the size of a planet grows, the 
corresponding feeding zone expands as well. Nevertheless, we 
could obtain the isolation mass of planet at a given orbital 
position by accumulating the total mass of the materials in 
the feeding zone of the original disk, 

/ M- \ 1 / 3 

M lso = 2na ■ 2F a ■ E p = 4na ■ C I -^-) a ■ E p (5) 
which gives, 



Mi ao oc a 2 • mV 3 ■ E„ a h 



3/4 



(6) 
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Figure 1. Snapshots of a run in the simulation of Group 1. The panels show the eccentricity versus semi-major axis for each surviving 
body at t = 0, 0.1, 1.0, 10, 50, and 400 Myr. The radii of each object are related to their masses, with radius <x m 1 / 3 . Two giants are, 
respectively, at 2.3 and 4.6 AU. Three terrestrial planets have formed at ~ 10 Myr and remain stable over ~ 400 Myr. 



Table 1. Orbital parameters of two gas-giants iGaudi et al.l 
l200Sf) . 



Planet 


a (AU) 


e 


Mass (M Jup ) 


OGLE-2006-BLG-109L b 


2.3 


0.001 


0.71 


OGLE-2006-BLG-109L c 


4.6 


0.11 


0.27 



Table 2. Simulation details of each group. 



Number 


Embryos 


Planetesimals 


Total Mass (M ffi ) 


Group 1 


30 


470 


10 


Group 2 


30 


470 


10 


Group 3 


50 


950 


15 


Group 4 


30 


470 


10 


Group 5 


30 


470 


10 



Therefore, we set the masses of planetary embryos and plan- 
etesimals following the relationship M em b r yo <x a 3 ^ 4 , pro- 
portional to the tota l mass in the embryo's feeding zone 
l|Ravmond et al.ll2004h . On the other hand, all the planetes- 
imals in our simulations were set with an equal mass of 0.017 
M®. 

The initial orbital elements of each embryo and plan- 
etesimal were randomly generated: argument of pericentre, 
longitude of the ascending node, and mean anomaly of each 
small object were randomly set between 0° to 360°; the ec- 
centricities range from to 0.02, while the inclination vary 
from 0° to 1°. For giant planets, to investigate terrestrial 
planetary formation in an inclined configuration, the incli- 
nation of the outer giant planet changes for each group but 
that of the inner gas-giant remains, in the meantime all other 
orbital data keep unchanged in the each initial run. 



As mentioned, one of the main goals of this work is to 
explore terrestrial planet formation in the late stage under 
the circumstance of highly-inclined planetary system, where 
the mutual inclination of the outer giant planet is taken into 
account. Consequently, we carried out five groups of simu- 
lations, containing 46 runs in total. The major difference 
between each group is the distribution range of the initial 
bodies. In all groups, two giant planets were set to emulate 
the OGLE-06-109L system in each run, with initial orbital 
parameters (M P , a, e p ) = (0.71 M iuv 2.3 A U, 0.001) and 
(0.27 Mjup, 4.6 AU, 0.11) (|Gaudi et all 120081 ). as shown in 
Table [J We adopt the stellar mass as 0.5 .M<n (|Gaudi et al.l 
2008), and assume its radius to be ~ 0.8 _Rq. For each group, 
the simply difference between individual runs is the inclina- 
tion of the outer giant, OGLE-2006-BLG-109Lc. Some de- 
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tails of each group are listed in Tabled which is also sum- 
marized as follows, 

• Group 1-5 runs. In each run, embryos and planetes- 
imals were distributed in the range 0.1 AU < a < 5.2 AU. 
The inclination of OGLE-2006-BLG-109Lc in each run is 
adopted from 0° to 40°, in increments of 10°. 

• Group 2-14 runs. In each run, embryos and planetesi- 
mals were in the range 0.3 AU < a < 5.2 AU. The inclination 
of OGLE-2006-BLG-109Lc in each run is from 0° to 10° in 
increments of 1°, and from 10° to 40° in increments of 10°. 

• Group 3-17 runs. In each run, embryos and planetesi- 
mals were in the range 0.3 AU < a < 5.2 AU. The inclination 
of OGLE-2006-BLG-109Lc in each run is from 0° to 10° in 
increments of 1°, and from 10° to 40° in increments of 5°. 

• Group 4-5 runs. In each run, embryos and planetesi- 
mals were distributed in the range 0.1 AU < a < 10 AU. he 
inclination of the outer giants in each run is from 0° to 40° , 
in increments of 10°. 

• Group 5-5 runs. In each run, embryos and planetes- 
imals were distributed in the range 0.3 AU < a < 10 AU. 
The inclination of the outer giants in each run is from 0° to 
40° , in increments of 10° . 

Both embryos and planetesimals were considered to 
gravitationally interact with each other in our simulations. 
Objects were allowed to collide, and we assume that they 
merge into a single body with no fragmentation after a col- 
lision. All 46 runs were performed using a hybrid symplectic 



algorithm provid ed by the MERCURY integration package 
i|Chambersll 19991) . Each run evolved for 400 Myr, with a time 
step length no more than 3.0 days (~ a twentieth of a pe- 
riod for the body at 0.3 AU). The Bulirsch-Stoer tolerance 
for all runs is of 10~ 12 . In most of the simulations, energy 
is conserved to better than 1 part in 10 -3 , and the angular 
momenta conserved in 10 -11 . Next, we will briefly introduce 
simulation outcomes. 



3 RESULTS 

From the simulations, we can note that a typical accre- 
tion scenario occurs in the l ate stage planet formation 
l|Chambersll200ll ; I Ji et~ai1l201ll '). In the beginning, the plan- 
etary embryos and planetesimals are quickly excited to 
highly-eccentric orbits due to dramatic perturbations arising 
from OGLE-2006-BLG-109Lb,c. Subsequently, frequent or- 
bital crossings emerge when the objects approach each other, 
which may result in violate collisions amongst planetesimals 
and embryos. Herein a large portion of the runs show that 
the planetesimal disk becomes quite turbulent within 1 Myr, 
while in some cases their chaotic period of movement could 
last a little longer, and most of the initial objects were re- 
moved by ejections or collisions. 

Fig. [T] illustrates that six snapshots of one run in Group 
1. The figure shows that three terrestrial planets emerged at 
~ 10 Myr and their orbits kept stable over subsequent evo- 
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Figure 3. Trace of the destination of all objects from 4 runs 
in Group 1. Each panel presents original versus final semi-major 
axes for all objects - except those scattered out of the system 
(a > 100 AU) - of which the mutual inclination for the giants are, 
respectively, 0°, 10°, 30°, and 40°. The vertical axis indicates the 
original semi-major axis of each body, while the horizontal axis 
indicates the final value (CS denote the central star). 



lution. Moreover, at the end of this run, there are also two 
survivals locked in the 1:1 Mean Motion Resonance (MMR) 
with OGLE-2006-BLG-109Lb, which reminds us of the Tro- 
jan/Greek population of asteroids at triangular points for 
Jupiter. In our simulations, survived Trojans are not rare to 
find, for they are simply unaccreted planetesimals. 

The dynamical evolution process within 10 Myr for the 
bodies plays a key role in creating members at final config- 
urations. As shown in Fig. O it represents the mass, semi- 
major axis, and eccentricity versus time for all final bodies. 
We observe that three terrestrial planets had finally made 
to reach about Earth-size mass in the stage. The eccentrici- 
ties of each body change violently during this chaotic stage, 
then they were damped down to ^ 0.1. In the final, the or- 
bits of three terrestrial planets remained stable afterwards. 
However, two planetesimals were trapped into 1:1 MMR 
with OGLE-2006-BLG-109b, and their orbital parameters 
are nearly unchanged over the timescale of 400 Myr. 

To summarize all outcomes, we note that the terrestrial 
planet formation is very common in the simulation, as a 
great number of runs may remain two or more planets even- 
tually. In addition, over 2/3 of the runs have 1-3 survival 
planetesimals at the 1:1 MMR with respect to OGLE-2006- 
BLG-109Lb. In the cases of the initial embryos and plan- 
etesimals extending 10 AU, the planetary accretion in the 
outer disk may still work. In such cases, a significant portion 
of objects could stay at the region 7-15 AU in the end. 



Figure 4. Trace of the destination of the objects from 4 runs in 
Group 4 — layout is the same as for Figure|3] In these runs, there 
are also some planetesimals survived in the outer planetesimals 
disk. 



3.1 Dependence on the outer highly-inclined giant 

With an increasing mutual inclination between two giants, a 
great many of smaller planetesimals either scattered out of 
the system, or had been directly engulfed by the central star. 
Therefore a few number of objects were left in the planetary 
disk, leading to the difficult of making terrestrial planets. 
Table [3] presents the statistic results over the simulations of 
Group 3, in which each run has an initial population of 1000 
bodies. The table summarizes the number of scattered bod- 
ies, those hit into giants and host star, the remaining mass 
and the accretion rate. For a moderate mutual inclination 
(im) of ~ 10° to 25°, the population of ejected bodies could 
be much more than that of i m < 5°. For a higher i m (> 25°), 
the number of ejected objects decreases. In this sense, a large 
relative inclination of two giants seems to be less effective in 
throwing out the embryos or planetesimals. However, on the 
other hand, in the highly-inclined configurations the pop- 
ulation that hit into the central star could be significantly 
enhanced. As shown in Table [3] with i m < 10° , simply a 
small number of planetesimals, say 20 - 40 objects, entered 
into the central star; while for a larger mutual inclination, 
the number could rise nearly linearly. Moreover, in the case 
of im = 30° , 40° , the bodies that originally lie in the outer 
disk might run into the central star over shorter timescale. 
In a word, the relatively inclined orbital configuration do 
play a significant role in the formation of terrestrial planets, 
inferring that the lower mutual inclination of two giants, the 
more effective for planetary accretion. 

Fig. [3] shows the destination of the initial objects in 
simulation of Group 1, except those were driven out of the 
systems. The vertical axis indicates the original semi-major 
axis of each body, while the horizontal line means the final. 
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Table 3. Statistic results for the simulations Group 3. The columns are, respectively, the mutual inclination, the number of ejected 
bodies, the objects collided with giant planets and central star, the remaining mass (in Earth-mass) and the surviving mass fraction for 
each run. 



Inclination (°) 


Ejection 


Giants 


Star 


Kemainmg Mass (Mq) 


Surviving Mass Fraction 





626 


69 


20 


3.404 


23% 


1 


669 


43 


26 


3.360 


22% 


2 


665 


54 


23 


3.450 


23% 


3 


664 


48 


20 


3.122 


21% 


1 


668 


37 


21 


3.526 


23% 


5 


686 


46 


28 


2.954 


20% 


6 


705 


49 


30 


2.701 


18% 


7 


707 


.39 


28 


3.230 


21% 


8 


725 


40 


39 


2.121 


14% 


9 


690 


53 


37 


2.460 


16% 


10 


748 


43 


31 


1.954 


13% 


15 


743 


43 


53 


1.491 


10% 


20 


746 


47 


76 


0.6314 


4% 


25 


752 


50 


89 


0.6184 


4% 


30 


711 


57 


107 


0.1441 


1% 


35 


699 


61 


105 


0.02275 


0% 


10 


685 


63 


120 


0.01137 


0% 



The figure shows that along with the augment of the mutual 
inclination of two jovian planets, the number of bodies that 
hit into the central star may increase notably. Especially, 
one may notice that the situation is more apparent in the 
inner planetesimal disk at 0.1 - 1.5 AU. Furthermore, Fig. [3] 
shows that in the case of i m = 30° , 40° , a large proportion 
of the objects in the inner disk could finally rush into the 
host star, leading to few residuals in the disk for further 
planetary accretion. This trend is in a good agreement with 
those as shown in Table [3] for Group 3. Again, Fig. [4] shows 
the destination of initial objects in Group 4, which reveals 
a similar trend like Fig. [3] A large number of bodies that 
finally run into the central star, may cause the deficiency in 
mass for the inner disk to produce more terrestrial planets, 
and a scarcity of survivals in the outer disk. 

As mentioned previously, a great many of bodies had 
been either driven out of the system, or collided with the 
giant planets or central star for inclined configuration, lead- 
ing to the difficulty in planetary accretion. Moreover, Fig. [5] 
shows the remaining bodies for each run from Group 3 over 
400 Myr, except for i m — 35°, 40° where no terrestrial planet 
formed. The figure clearly reveals that the mass of created 
planets decreases as the mutual inclination grows. As shown 
in Fig. O the i m = 30° run illustrates that a terrestrial 
planet with mass of ~ 0.1M® was finally formed at ~ 0.8 
AU, indicative of inefficiency for planetary accretion. To sum 
up, simulation results demonstrate that terrestrial planetary 
formation is extremely related to the mutual inclined con- 
figuration of the giant planets in the system. 

3.2 Dependence on the orbital distribution of 
initial bodies 

The different distribution ranges of initial bodies also play a 
major part in yielding final members in the system for each 
simulation. Table [4] introduces five contrast groups, and in 
each run the number of the initial bodies is 500, with a total 
mass of 10 M®, but they are distributed in various initial 
orbits. As shown in Table(4] in each group - in which all runs 



have equal mutual inclination relative to giant planets - runs 
with an initial disk of an inner edge at 0.1 AU seem to have 
advantage of occupying more survivals of terrestrial planets 
and planetesimals, with a larger total mass of remaining 
bodies and efficient accretion rate of the terrestrial planets. 

Fig. [S]shows the outcomes of ten runs from Group 1 and 
2, to compare the overall results for different initial distribu- 
tion ranges. Similar to Table|4j both the number and mass of 
formed terrestrial planets in Group 1 (0.1 AU < a < 5.2 AU) 
is larger than those in Group 2 (0.3 AU < a < 5.2 AU). Con- 
sidering the smaller semi-major axis of two giants in this 
system, planetary embryos and planetesimals in the outer 
region of the planetesimals disk would suffer severe gravita- 
tional influence from two giant planets; hence, the planetes- 
imals, locating at 1 ~ 7 AU could not remain stable under 
dramatic perturbation from them. The only exception is the 
1:1 MMR region with respect to the OGLE-2006-BLG-109b, 
while no such Trojan was found for the outer giant planet. 
This may explain an almost empty disk for all final configu- 
rations, ranging from 1 ~ 7 AU, in our simulations. On the 
contrary, the inner part of the disk is strongly favor of pro- 
ducing terrestrial planets, because of the relatively higher 
density of the population and stronger influence of the cen- 
tral star. In this sense, these cause the inner region of the 
planetary disk, from 0.1 to 0.3 AU, to be very important for 
planetary accretion for terrestrial planets. 

The above feature is also reported in much wider distri- 
bution cases, where the embryos and planetesimals originally 
located at an outer edge of disk, extending to 10 AU. Figure 
[7] shows the terrestrial planetary accretion becomes more 
difficult when the mutual inclination of the gas-giants in- 
creases. One may note that in the final evolution there exist 
several survivals in the outer region of the disk up to 7 AU. 
In Figure [8] the runs of 0.1 ~ 10 AU configuration would 
relatively have less and smaller final planets than those in 
the 0.3 ~ 10 AU cases. 

In addition, in the cases of initial distribution up to 10 
AU, there also exist several survivals in the outer realm of 
the planetary disk at 7 AU. Some survivals experienced sev- 
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Figure 5. All the final configurations of the runs that have terrestrial planets formed in Group 3. The panels show the eccentricity 
versus semi-major axis for each surviving body at the end of each run. The radii and colors of each object are related to their masses, 
with radius <x m 1 / 3 . The mutual inclination of the two giants is labeled in each panel. Note that the mass of formed planets decreases 
as the mutual inclination between the giants increase. 



eral accretion, at a limited rate, in comparison with those in 
the inner disk; others, like Mar-sized protoplanets, emerged 
around the region of 7.5 ~ 11 AU in the end of simulations. 
Compared with our solar system, these formed bodies may 
be in a rescaled position of Uranus, Neptune and the Kuiper 
belt objects, implying that Neptunian-mass planet may be 
formed if more bodies were initial placed there and abundant 
mass were considered over much longer timescale evolution. 

3.3 Planetary formation in the habitable zone 

As mentioned, the OGLE-2006-BLG-109L was the first dou- 
ble plane t system, discovere d by gravitational microlensing 
method (|Gaudi et al.l l2008h . Two companions occupy the 
masses of ~ 0.71 and ~ 0.27Mj, respectively, and each or- 
bits its central star at about 2.3 and 4.6 AU. The mass of 
the star is ~0.5 Mq, and this system is resemblance to a 
rescaled solar system, for the mass ratio and the separation 
between them is comparable to those of Jupiter and Saturn. 
The habitable zone around a star is defined as the region 



where a terrestrial planet with N2-CO2-H2 O atmosphere 
could maintain liquid oceans on its surface jKasting et al.l 
Il995t ). As for this system, the HZ is a bout 0.25~0.36 AU 
from the star (|Migaszewski et al.l [20091 ) . Since the OGLE- 
2006-BLG-109L planetary system bears very close similarity 
to solar system, one may wonder to know whether it is pos- 
sible to harbor other planets, es pecially for potential Earth - 
like planets in the HZ. Recently. iMalhotra fc Mintonl (|200St ) 
analyzed secular dynamics for OGLE-2006-BLG-109L and 
presented a stable model to explain the e xistence of two sup- 
posed terre strial planet s in th e system. iMigaszewski et ahl 
l|2009h and IWang et all l|2009l ') further showed that Earth- 
size bodies might be formed in the HZ of the system. 

From our simulations, we note that more than 40% (19 
out of 46) of the total runs have formed planet s in the HZ, 
within 0.25 ~ 0.36 AU l|Migaszewski et al.|[2009l '). In fact, we 
find that in certain configurations with a smaller mutual i 
(^ 10°) may have a great probability, e.g., 16 out of 22 runs 
formed a terrestrial planet in the HZ. Fig. [9] shows all final 
configurations of 19 runs have produced habitable planets 
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Table 4. Statistics over 20 selected runs from four different groups, considering the mutual inclinations. The name of each run denotes 
the initial distribution range and the relative inclination between two giants, e.g., 10-0.1-00 indicates a run for embryos and planetesimals 
initially ranging from 0.1 to ~ 10 AU, with a 0° mutual inclination for two giant planets. In the following columns, the table lists the 
final number of remaining terrestrial planets and planetesimals (referred as TP and PL, respectively), the total remaining mass (in 
Earth-mass) and the surviving mass fraction of formed terrestrial planets. 



Simulation 


Initial 


Final(TP + PL) 


Remaining Mass (Mg) 


Surviving Mass Fraction (TP) 


i n n 1 nn 
1U-U. 1-UU 


OUU 




q non 

o.uyu 


1 i /o 


10-0.3-00 


500 


2+13 


2.832 


27% 


5.2-0.1-00 


500 


3+2 


2.852 


28% 


5.2-0.3-00 


500 


2+2 


1.900 


19% 


1 n n 1 in 
1U-U.1-1U 


OUU 


Till) 


Z.DoU 


24 /o 


10-0.3-10 


500 


2+10 


1.585 


13% 


5.2-0.1-10 


500 


3+2 


2.469 


24% 


5.2-0.3-10 


500 


2+2 


1.604 


16% 


10-0.1-20 


500 


3+9 


1.946 


18% 


10-0.3-20 


500 


1+6 


0.777 


7% 


5.2-0.1-20 


500 


2+2 


1.864 


18% 


5.2-0.3-20 


500 


3+1 


0.931 


9% 


10-0.1-30 


500 


3+8 


1.174 


11% 


10-0.3-30 


500 


0+8 


0.381 


3% 


5.2-0.1-30 


500 


2+1 


1.334 


13% 


5.2-0.3-30 


500 


1+1 


0.210 


2% 


10-0.1-40 


500 


2+4 


0.436 


4% 


10-0.3-40 


500 


0+6 


0.222 


1% 


5.2-0.1-40 


500 


2+0 


1.363 


14% 


5.2-0.3-40 


500 


1+3 


0.229 


3% 



at the end of numerical simulation, where Run 1 ~ 8 are 
adopted from Group 2, Run 9 ~ 16 from Group 3, and Run 
17, 18 and 19 from Group 1, 4 and 5, respectively. Obviously, 
as illustrated in Fig. [9] most of the runs are from Group 2 
and 3, which have an initial distribution for embryos and 
planetesimals ranging from 0.3 - 5.2 AU. 

The common feature in these runs is indicative of the 
final configuration to consist of two terrestrial planets, with 
one in the HZ and the other in outer region from 0.4 to 
1.0 AU. Typically, the planet in the HZ seems to be more 
massive (0.75 ~ 2.6M©), while the outer body has a smaller 
mass (0.27 ~ 1.81M®). In addition, there also exists some 
configurations (totally 6 runs) that three terrestrial planets 
formed. In these cases, the other two planets do not reside in 
habitable zone. However, no configuration of a single habit- 
able planet has been found in the final simulations, although 
it was considered to be a stable model under the combina- 
tions of a low eccentric or bit and some proper orb ital pa- 
rameters of jovian planets (jMigaszewski et aDl2009h . 

The uniform characteristics of the formation of hab- 
itable terrestrial planets in our simulations over 400 Myr, 
could be served as a significant evidence of the stability of 
the HZ in the OGLE-2006-BLG-109L system. 



4 DISCUSSION 

The OGLE-2006-BLG-109L system has drawn many re- 
searche r's attention as it is resembling to a rescaled solar 
system. IWang et alJ (2009) showed that the stability of the 
region is within a ^ 1.5 AU and a ^ 9.7 in their numer- 
ical simulations, and also pointed out the HZ is broader 



enough for a stable Earth- mass planet. From the runs, we 
find that terrestrial planets could emerge in the regions 
where 0.1 < a < 1.13 AU or a > 7.5 AU. For 1.13 < a < 1.5 
AU, there may remain stable for an additional planet as sug- 
gested bv lWang et all (|2009l ). where a great many of initial 
bodies are ejected due to strong stirring of two jovian plan- 
ets, and they either scatter out of the system or directly run 
into the central star. As a result, simply a few number of 
objects were left as accretion stuff for further planetary for- 
mation, consequently none of the planets has finally formed 
in this region. Herein we also confirm the stability of the HZ 
in this system: 19 out of total 46 runs have a final planet 
locating in 0.25 to 0.36 AU, which is in good agreement with 
theirs. 

iMalhotra fc Mintonl (2008) found that secular reso- 
nance arising from the massive planets in the system could 
drive an Earth-mass planet out of the HZ, and an addi- 
tional inner planet, with mass larger than 0.3M® at a ^ 0.1 
AU, is required to maintain the habitability of the OGLE- 
2006-BLG-109L system. Therefore, they suggest a possible 
configuration with at least two additional planets in order 
to support a potentially habitable Earth-like planet in this 
system. However, in the 19 simulations that formed a habit- 
able planet, we show that most of final configurations bear 
one larger planet (0.75 ~ 2.6Mq) in the HZ, and an addi- 
tional one or two less massive (0.27 ~ 1.81M©) planets in 
the outer region ranging from 0.4 to 1.0 AU. Simply one 
run in our work may have a similarity to the configuration 
as mentioned by IMalhotra fc Mintonl ()2008f ) . where three 
planets with parameters (M P , a) = (0.1 M®, 0.0385 AU), 
(2.05 M e , 0.3044 AU), and (0.55 M©, 0.7951 AU). Our nu- 
merical results a bit differ from their theoretical analysis 
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Figure 6. The comparison of final configurations of 5 runs in Group 1 (where 0.1 AU < a < 5.2 AU for initial objects) and Group 2 
(0.3 AU < a < 5.2 AU), and the layout is the same as FigureO Note that in the case of Group 1 more terrestrial planets could form. 



in the initial distribution of inner boundary (where we set 
to be 0.1 AU) for planetesimals and embryos. However, from 
viewpoint of habitability, our simulations are in consistent 
with their suggestion that at least two additional planets are 
required to support a potentially habitable Earth-like planet 
in this system. 

In the simulations, the configuration s with a single sta- 
ble h abitable planet as mentioned by iMigaszewski et all 
(2009) do not occur, where a terrestrial planet could sur- 
vive to move on low-eccentric orbit for the combination of 
the apsidal angle of jovian orbits, eccentricities and semi- 
major axes consistent with the observations. However, our 
work differs from those in that in our simulations we adopt 
the highly-inclined configuration between two jovian plan- 
ets, and we also have set a much larger total mass an d more 
initial bodies than those of lMigaszewski et alj |2009l ) , where 
they simply considered less bodies - 50 Moon-sized planetary 
embryos between 0.2 - 2 AU, leading to a great difficulty in 
producing more terrestrial planets. Nonetheless, our work 
is in concurrence with their study in the fact that most of 
the runs formed a habitable planet have an additional body 
with a similar mass in the region 0.6 ~ 0.8 AU. 



Recently, IWang fc Zhoul (|2010h investigated the forma- 
tion of terrestrial planets in the early stage of planet for- 
mation at the time the gaseous disk has not dissipated yet. 
They indicated that the formation and the final mass of the 
habitable planet in the OGLE-06-10L system depends on 
the stellar accretion rate M, which means the rate at which 
the star accretes gas from the disk, and the speed of Type 
I migration, which is more significant. In this work, we take 
into account formation of terrestrial planets in the late stage 
after the gas disk has dissipated. We concentrate on the var- 
ious distribution range of the initial bodies and mutual incli- 
nation between two gas giants. Generally, the higher relative 
inclination, the harder the planetary accretion. This is be- 
cause the relative highly-inclined outer planet may disturb 
the embryos and planetesimals in a more effective manner, 
which results in a large number of objects thrown out of 
the system or run into the central star, and then a limited 
number of bodies are left in the disk for further planetary 
accretion. 

With a relatively higher inclination, a large amount of 
the bodies would be excited to be swallowed by the central 
star, especially those in the inner region 0.1 to 1.5 AU. More- 
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Figure 7. Final configurations of 5 runs in Group 4. The layout is the same as Figure[5l but the horizontal line is in a logarithmic scale. 
The figure shows the terrestrial planetary accretion becomes difficult in the case of the increasing mutual inclination of the gas-giants. 
In the final, there exist several survivals in the outer region of the disk extending to 7 AU. 



over, the total mass of the scattered bodies also increases as 
the mutual inclination does. In a word, the above-mentioned 
two factors do no good to form terrestrial planets for the 
highly-inclined configurations. 

Given two giants in this system with close-in orbits, 
the embryos and planetesimals in the inner disk were un- 
der strong gravitational influence arising from the jovian 
planets. Our results show that the total accretion mass in 
the runs of the initial distribution at 0.1 AU is much larger 
than those distributed at 0.3 AU. In this sense, it may imply 
that the location of the inner planetesimal disk, as well as 
the number of the initial bodies, plays an important part 
in the planetary accretion. On the other hand, embryos and 
planetesimals between 1 ~ 7 AU could not survive so long 
owing to dramatic stirring of two gas-giants, except those at 
1:1 MMR location with the inner planet, which reminds us 
Trojan/Greek Group of Jupiter. In two groups, the accretion 
does still work in the outer region beyond 7 AU, with a mod- 
erate rate and limited speed. Figs. [7] and [5] show that several 
Mar-sized protoplanets emerge between 7.5 to 11 AU over 
400 Myr. In addition, there are also many survival planetes- 
imals at 7 ~ 16 AU. Suppose a larger number and a longer 
integration time were chosen at the beginning of these runs, 
it may assume that Neptune-like planets might appear in 



the end due to further planetary accretion in the outer disk, 
similar to the outer structure of our solar system. 



5 SUMMARY 

In this work, we have carried out several groups of simula- 
tions, with 46 runs in total, to investigate the planet forma- 
tion in the late stage for the OGLE-2006-BLG-109L system. 
We may summarize the main results as follows: 

Firstly, the simulation outcomes show that it is quite 
common for terrestrial planet or possible habitable plan- 
ets to finally form residing in 0.25 ~ 0.36 AU. In addition, 
our work may imply that the OGLE-2006-BLG-109L sys- 
tem bears a great probability of harboring one or additional 
terrestrial planets. 

Secondly, the comparison of results for different groups 
suggests that the efficiency of planetary accretion may de- 
crease as the relative inclination of two giant planets in- 
creases. 

Finally, our results further show that embryos and plan- 
etesimals in the 0.1 - 0.3 AU may play an important role in 
planetary accretion. Moreover, we find that the bodies ini- 
tially locate in the ~ 1.1 AU - 7 AU of planetary disk are 
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Figure 8. All the final configurations of 5 runs in Group 4, the layout is the same as Figure In Comparison with Figure [7] the runs 
of 0.1 ~ 10 AU configuration would relatively have less and smaller final planets than in the 0.3 ~ 10 AU distribution cases. 



mostly unstable over the dynamical evolution, which leads 
to no survivals for embryo or planetesimal at this range in 
all simulations, except those at ~ 2.2 AU close to 1:1 MMR 
with OGLE-2006-BLG-109Lb. On the other hand, the plan- 
etary accretion could continue to happen but keep a much 
slower pace in the outer planetary disk for those runs with 
initial semi-major axes extending 10 AU, where a significant 
number of residuals eventually remain in the system in the 
simulations. 

Currently, the Kepler mission is surveying 156,000 faint 
stars for transiting p lanets as small as Earth, e.g., Kepler- 9 
ijHolman et all l2010h and Kepler- 11 (jLissauer et al.ll201ll ). 
Our results here imply the great probability to occupy sta- 
ble terrestrial planets in the OGLE-2006-BLG-109L system, 
which should be carefully examined by forthcoming obser- 
vations. With the help of high precision of ground-based 
and spaceborne measurements on the exoplanets, we might 
expect the discovery of solar system analogues in the near 
future. 
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